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ABSTRACT 

We report new continuum observations of fourteen z~6 quasars at 250 
GHz and fourteen quasars at 1.4 GHz. We summarize all recent millimeter 
and radio observations of the sample of the thirty-three quasars known with 
5.71 < z < 6.43, and present a study of the rest frame far- infrared (FIR) proper- 
ties of this sample. These quasars were observed with the Max Plank Millimeter 
Bolometer Array (MAMBO) at 250 GHz with mJy sensitivity, and 30% of them 
were detected. We also recover the average 250 GHz flux density of the MAMBO 
undetected sources at 4cr, by stacking the on-source measurements. The derived 
mean radio-to-UV spectral energy distributions (SEDs) of the full sample and 
the 250 GHz non-detections show no significant difference from that of lower- 
redshift optical quasars. Obvious FIR excesses are seen in the individual SEDs 
of the strong 250 GHz detections, with FIR-to-radio emission ratios consistent 
with that of typical star forming galaxies. Most 250 GHz-detected sources follow 
the LpiR-Lboi relationship derived from a sample of local IR luminous quasars 
(Ltr > 1O 12 L ), while the average LprR/Lboi ratio of the non-detections is con- 
sistent with that of the optically-selected PG quasars. The MAMBO detections 
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also tend to have weaker Lyct emission than the non-detected sources. We discuss 
possible FIR dust heating sources, and critically assess the possibility of active 
star formation in the host galaxies of the z~6 quasars. The average star forma- 
tion rate of the MAMBO non-detections is likely to be less than a few hundred 
M Q yr _1 , but in the strong detections, the host galaxy star formation is probably 
at a rate of > lO 3 M yr _1 , which dominates the FIR dust heating. 

Subject headings: galaxies: quasars - galaxies: active - submillimeter - galaxies: 
starburst - galaxies: high-redshift 

1. Introduction 

More than thirty quasars at z~6 have been discovered (e.g., Fan et al. 2000, 2001a, 
2003, 2004, 2006a; Jiang et al. 2007a; Willott et al. 2007). These objects are the earliest 
massive black hole accretion systems known (Jiang et al. 2006; 2007b), seen at an epoch 
close to the end of cosmic reionization (Fan et al. 2006b). They provide us with unique 
information on both the growth of supermassive black holes (SMBHs) and the formation of 
massive galaxies when the age of the universe was < 1 Gyr. 

There are fundamental relationships between SMBH mass and bulge stellar mass/velocity 
dispersion in the local universe (e.g., Tremaine et al. 2002; Marconi & Hunt 2003; Hopkins 
et al. 2007), indicating that the formation of SMBHs and their spheroidal hosts are coupled. 
Active galactic nuclei (AGNs) buried in dusty starburst environments have been discovered 
in samples of low redshift ultraluminous infrared galaxies (ULIRGs), which are believed to 
be a transition phase between starburst mergers and typical optically bright AGNs (Sanders 
et al. 1988; Wu et al. 1998; Zheng et al. 2002). Studies of large samples of galaxies and 
AGNs also provide clues about the 'downsizing' process in galaxy and SMBH formation, 
whereby the formation of the most massive systems occur at early epochs (Cowie et al. 
1996; Heckman et al. 2004; Kauffmann & Heckman 2005). 

Examples of massive quasars hosted by interacting systems or bright submillimeter 
galaxies have been found at redshifts greater than 4, such as BR 1202-0725 at z=4.7 and 
BRI1335-0417 at z=4.4 (Omont et al. 1996a, 1996b; Benford et al. 1999; Beelen et al. 2006). 
These objects are all characterized by strong molecular CO line emission (Guilloteau et al. 
1997, 1999; Carilli et al. 2002; Solomon & Vanden Bout 2005; Riechers et al. 2006) and FIR 
(Guilloteau et al. 1999) and radio (Momjian et al. 2005, 2007) emission originating from a 
starburst, with implied star formation rates of a few thousand M Q yr _1 . The results suggest 
that the galaxies are at an earlier evolution stage than are typical optically bright quasars 
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in which the central AGN dominates the emission from X-ray to radio. 

These studies have been extended to z~6 quasar sample with sensitive submillimeter, 
millimeter, and radio telescopes (Priddey et al. 2003b; Bertoldi et al. 2003; Petric et al. 
2003; Carilli et al. 2004; Wang et al. 2007). About 1/3 of these z~6 quasars were detected 
at millimeter wavelengths, at mJy sensitivity (Priddey et al. 2003b; Wang et al. 2007). The 
millimeter detections imply FIR luminosities of 10 12 ~ 10 13 L Q and dust masses of > 10 8 M Q 
in the quasar host galaxies (Bertoldi et al. 2003a, Beelen et al. 2006). Such huge dust 
masses require rapid metal and dust enrichment via a starburst in the early evolution of 
these galaxies (Bertoldi et al. 2003a; Walter et al. 2003; Venkatesan et al. 2006; Li et al. 
2007b; Dwek et al. 2007). 

The heating sources of the FIR-emitting dust in the quasar systems at z~6 have become 
a key question. The reprocessed emission from star formation-heated dust can provide a 
direct estimate of the star formation rate, thus constraining the bulge building stage in these 
quasar hosts. Observations at submillimeter wavelengths imply that the FIR emission in the 
two sources with the strongest millimeter detections, SDSS J114816. 64+525150. 3 at z=6.42 
(hereafter J1148+5251) and SDSS J092721.82+200123.7 at z=5.77 (hereafter J0927+2001), 
is from dust components with temperatures of 50 to 60 K (Beelen et al. 2006; Wang et 
al. 2008a). Large amounts of highly excited molecular CO were also detected in the host 
galaxies of these two sources (Bertoldi et al. 2003b; Walter et al. 2003; Carilli et al. 2007), 
as well as strong [C II] 158 /im ISM gas cooling line emission in J1148+5251 (Maiolino et al. 
2005). These results suggest that star formation plays an important role in the heating of 
FIR-emitting warm dust. The implied star formation rate is > lO 3 M yr _1 , which argues 
for active bulge building in these two z~6 quasars. 

We have a long standing program to study the dust and gas emission from the z~6 
quasars (e.g., Wang et al. 2007, 2008a). In this paper, we report new observations and 
present an FIR emission study of the full sample of thirty-three z~6 quasars discovered to 
date, aiming (i) to give a general view of the FIR emission properties of the sample, and 
(ii) to further constrain the dust heating and star forming activity in the quasar hosts. We 
describe the full quasar sample, the new observations, and summarize the current millimeter 
and radio results in section 2. The analysis of the full sample is given in Sections 3, 4, and 5. 
We present a discussion of star formation in Section 6, and give a brief summary in Section 
7. We adopt a A-model cosmology with Ho = 71km s _1 Mpc -1 , Qm = 0.27 and = 0.73 
throughout this paper (Spergel et al. 2007). 
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2. Sample and observations 

2.1. The sample of z~6 quasars 

Thirty-three quasars have been discovered at z~6 (5.71 < z < 6.43, see Table 1 and 
2). Twenty- two of these objects were selected from ~8000 deg 2 of the Sloan Digital Sky 
Survey (SDSS, York et al. 2000) imaging data (Fan et al. 2000, 2001a, 2003, 2004, 2006a, 
2008, in prep.; Goto 2006). These objects typically have rest frame 1450A AB magnitudes 
m i450 < 20. The SDSS imaging survey is using a dedicated 2.5m telescope and a large 
format CCD camera (Gunn et al. 2006) at the Apache Point Observatory in New Mexico to 
obtain images in five broad bands (u, g, r, i and z, centered at 3551, 4686, 6166, 7480 and 
8932 A, respectively; Fukugita et al. 1996) of high Galactic latitude sky in the Northern 
Galactic Cap. The imaging data are processed by the astrometric pipeline (Pier et al. 2003) 
and photometric pipeline (Lupton et al. 2001), resulting in astrometric calibration errors of 
< 0.1" rms per coordinate, and photometric calibration to better than 0.03 mag (Hogg et 
al. 2001, Smith et al. 2002, Ivezic et al. 2004, Tucker et al. 2006). These data have been 
made available to the public in a series of data releases (EDR: Stoughton et al. 2002; DR1 
- DR6: Abazajian et al. 2003, 2004, 2005, Adelman-McCarthy et al. 2006, 2007, 2008). 

Five additional z~6 quasars were discovered in the SDSS Deep Equatorial Stripe ('Stripe 
82') by Jiang et al. (2007a). These five sources were selected from deep optical imaging of 260 
deg 2 . These objects are about one magnitude fainter in the optical (i.e., 20 < m 145 o < 21.5) 
than the primary SDSS sample. One source in this sample, J020332. 39+001229. 3 (hereafter 
J0203+0012), was also independently discovered in the UKIRT Infrared Deep Sky Survey 
(Venemans et al. 2007). 

Four new quasars at z > 6 have been published (Willott et al. 2007) from the Canada- 
France High-z Quasar Survey (CFHQS). These quasars were selected from imaging of 400 
deg 2 ; three of them have m 145 o > 21. The fourth source, CFHQS J1509-1749, is comparable 
in brightness to the faintest source in the primary SDSS sample, but the optical spectrum 
is much steeper than that of the typical quasars, indicating significant dust reddening. 

Another two z~6 quasars were discovered in the NOAO Deep Wide-Field Survey (Jan- 
nuzi et al. 2004) using optical, infrared, and radio data (McGreer et al. 2006; Cool et al. 
2006). One of them, J1427+3312 at z=6.12, was selected by McGreer et al. (2006) from 4 
deg 2 by matching infrared and radio sources to optical counterparts. It is the brightest radio 
source in the z~6 quasar sample. The UV to optical spectrum of this object is likely to be 
reddened by dust, and the rest frame 1450A magnitude derived from near- infrared data is 
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mi45o = 20.33j. The other quasar was infrared selected from Spitzer mid-infrared data, and 
identified in the AGN and Galaxy Evolution Survey (Cool et al. 2006) with m 14 5 = 20.62. 

2.2. New observations 

Twenty-two sources in the z~6 quasar sample have published 250 GHz continuum ob- 
servations with the Max-Planck Millimeter Bolometer Array (MAMBO, Kreysa et al. 1998) 
on the IRAM 30-m telescope (Bertoldi et al. 2003a; Petric et al. 2003; Wang et al. 2007; 
Willott et al. 2007; Priddey et al. 2008), and nineteen sources in the sample have published 
1.4 GHz continuum observations with the VLA (Petric et al. 2003; Carilli et al. 2004; Wang 
et al. 2007; McGreer et al. 2006). We report our new MAMBO and VLA observations 
of totally nineteen z~6 quasars in this paper. These sources are unobserved or marginally 
detected in the previous observations, and we present new data for them at either 250 GHz 
or 1.4 GHz. We describe the new observations below and list the measurements for the 
nineteen quasars in Table 1. We also summarise the previous observations for other z~6 
quasars in Table 2. 

We present new 250 GHz observations for fourteen z~6 quasars, including eleven sources 
that have never been previously observed. The new observations have been obtained with the 
MAMBO-II 117-element array in Winter 2006-2007, we adopted the same observing mode as 
was used in the previous observations (Wang et al. 2007), i.e., doing chopping photometry 
at 2 Hz with 32" throw in azimuth. We typically spent 1-3 hours on each source to achieve 
the sub-mJy rms noise level reached in previous observations. The data were reduced with 
the standard MOPSIC pipeline (Zylka 1998). Three sources, J1044-0125, J1425+3254, and 
J2054-0005, were detected from the new data at > 4cr level. The MAMBO detection of 
J1044-0125 was recently reported in Priddey et al. (2008) from an independent observation 
with a 250 GHz flux density of 2.5 ± 0.6. We observed this source down to a lower rms noise 
level, and our result, 1.8 ± 0.4 mJy, is a bit lower compared to the value in Priddey et al. 
(2008) but consistent within the error limits (see Table 1). By including the data from the 
new observations, we also confirmed the previous 3a detection in the source J0818+1722 
(Wang et al. 2007) with a lower rms of 0.4 mJy. 

New 1.4 GHz data were obtained with the VLA for fourteen quasars, including thir- 
teen previously unobserved sources and a published marginally detected source, J0927+2001 
(3.2a, Wang et al. 2007). The new observations were carried out in 2006 and 2007 with the 
array in the A or B configurations, down to an rms noise level of ~ 20 /xJy for most of the 



Merived from the absolute magnitude of M1450 = —26.4 presented in McGreen et al. (2006) 
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sources. There is a ~1 Jy source in the field of the quasar J0303-0019, leading to a high rms 
of 66 //Jy in the image. All the data were reduced using AIPS. One source, J0203+0012, was 
detected by our new observation with a 1.4 GHz flux density of S u = 195 ± 22 //Jy, and the 
previous detection of J0927+2001 was confirmed at > 4<r level. 

There are fairly bright and extended radio sources in the 1.4 GHz field images of the 
CFHQS quasars J0033-0125 and J1641+3755 (Figure 1), though these two quasars are not 
detected. One double-lobe radio source is found about 9' away from the optical position 
of the quasar J0033-0125, with integrated 1.4 GHz flux density of ~ 23mJy. This radio 
source was also detected in the NVSS survey (Condon et al. 1998) and FIRST survey 
(White et al. 1997), and there is no published deep optical data for this area yet. Five 
radio sources with integrated flux densities > 10 mJy are seen in the field of the source 
J1641+3755, including extended double radio sources, within 10' from the optical quasar 
position. According to the FIRST 1.4 GHz source counts (White et al. 1997), about 1.6 
source with fi.4GHz > 10 mJy is expected within this amount of sky area. The brightest radio 
source in the field of J1641+3755 was identified as a radio-loud AGN with a spectroscopic 
redshift of z=0.162 (Best et al. 2005). Further optical imaging and spectrography will 
address the question of whether the bright radio sources in the fields of the two CFHQS 
quasars is due to foreground galaxy clusters. We notice that bright extended radio sources 
and overdensities were previously reported in the 1.4 GHz observations of the two SDSS 
z~6 quasars J1030+0524 and Jl 148+5251 (Petric et al. 2003; Carilli et al. 2004). Carilli 
et al. (2004) investigated the SDSS data of the J1148+5251 field, and however, there is no 
evidence for dense foreground clusters or gravitational lensing. 



2.3. Summary of the millimeter and radio results 

In summary, all the thirty-three source in the z~6 quasar sample have been observed by 
MAMBO at 250 GHz. The observations have rms values in the range of 0.4 — 1.1 mJy with 
a median value of about 0.6 mJy, and ten sources are detected (see Table 1 and 2). This 
yields a detection rate of 30 ± 10% which is consistent with the submillimeter and millimeter 
detection rates of optically bright quasars at redshifts 2 and 4 (Priddey et al. 2003a; Omont 
et al. 2001; 2003; Carilli et al. 2001). 

Thirty-two z~6 quasars have 1.4 GHz radio continuum observations by the VLA, with 
a median rms level of 20 //Jy, Ten out of the thirty-two sources have been detected, including 
four 250 GHz detections (see Table 1 and 2). Two of the radio detections, J0836+0054 and 
J1427+3312, have flux densities > lmJy. These two sources, together with the new radio 
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detection J0203+0012, have radio loudnesgj R = fi/,5GHz/fi/,44oo > 10, and hence are radio 
loud according to the definition by Kellermann et al. (1989). 

Data at other wavelengths for the analysis in this paper were collected from the liter- 
ature, including optical data from the original papers (Fan et al. 2000; 2001a; 2003; 2004; 
2006a; 2008, in prep; Goto 2006; McGreer et al. 2006; Willott et al. 2007; Jiang et al. 2008), 
near-infrared data from Jiang et al. (2006), submillimeter data from Priddey et al. (2003b), 
Robson et al. (2004) and Beleen et al. (2006), and the millimeter continuum data from the 
PdBI CO observations (Carilli et al. 2007; Wang et al. 2008b, in prep.). 



3. The UV-to-radio SEDs of the z~6 quasars 

Using the 250 GHz and 1.4 GHz continuum measurements, we derive the average FIR 
and radio emission of these z~6 quasars and compare it to the typical quasars at low red- 
shift. The calculations are performed for three groups: (I) the whole sample of thirty-three 
sources, (II) the subsample of ten sources detected with MAMBO at > 3cr, and (III) the 
subsample of twenty-three sources not detected with MAMBO. Using the MAMBO 250GHz 
data, we calculate the weighted average 250GHz flux densities for each group, namely: 

,r \ ^Wif250GHz,i 
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where f250GHz,i is the MAMBO observed flux density for each source, and 02^GHz,i is the la 
rms. For the calculation of (f250GHz) of the whole sample and the subsample of 250 GHz 
non-detections, we weighed the data by l/ofsoGHz.i (equation (2)). However, equal weights 
were adopted (i.e. Wi = 1) when calculate (f 2 50GHz) for the 250 GHz detections, as the 
observed rms values spread in a smaller range (0.36~0.75 mJy) and we do not want to bias 
the average toward the weakest sources. The average flux density ratios between 250GHz 
and the rest frame 1450A are then calculated as (f 2 50GHz)/(fi45o), where the average 1450A 
flux density (fi45o) is derived from the AB magnitudes presented in the discovering papers 



2 The 4400A and 5 GHz flux densities are derived from the 1450A magnitude and the observed 1.4 GHz 
flux density, assuming power law spectra indices of -0.5 in the optical and -0.75 in the radio, respectively. 
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(see Section 2.1). The average radio emission and 1.4GHz-to-1450A ratios for each group 
are calculated using the same method. 

The average FIR luminosities (42.5/xm - 122.5/xm, rest frame) for each group can be 
estimated^] as < L FIR >= 2.34 x 10 12 (< f250GHz > / m Jy)L , assuming a dust temperature 
of 47 K and an emissivity index (3 of 1.6 (Beelen et al. 2006). Based on this estimation, the 
mean FIR-to-radio ratios for each group can be calculated as (Helou et al. 1985): 

q = log ( (Lfir) V log ( (Ll4GHz) ) (4) 

The rest-frame 1.4 GHz radio luminosity densities (L L4 gh z ) are derived with the average 
radio flux densities, assuming a power- law radio spectrum of i v ~ //- - 75 (Condon 1992). 
The q values can be increased by ~0.2 if we assume a flater radio spectrum of i v ~ z/~ 05 . 
The mean flux densities, luminosities, and q values are listed in Table 3. We exclude the 
three radio loud sources in the first and the third groups, and present the corresponding 
results in the second row of each group in Table 3. The ten MAMBO detected sources in 
the second group are all radio quiet (according to the definition in Kellermann et al. 1989). 



3.1. The mean SED 

The mean 250 GHz flux density is 1.26 ± 0.10 mJy for the whole sample, and 2.30 ± 
0.15 mJy for the MAMBO detections. The mean 250 GHz flux density of MAMBO non- 
detections is recovered by averaging all the on-source pixel values. A ~ 4<r signal of 
0.52 ± 0.13 mJy is obtained for this group, which is about one quarter of the value of the 
typical MAMBO detections. We ran the stacking with the non-detection group several times, 
omitting one source each stack, to see if this 4<r signal is due to a single source that is just 
below the detectability. The test shows stacking flux densities in the range of 0.47 to 0.56 
mJy with rms values of ~0.13 mJy. Thus it is unlikely that the average flux density is 
dominated by one object. No significant differences are found in the (f250GHz) values when 
the radio- loud sources are excluded in each group. 

The average 1.4 GHz observed flux densities are 27±4to37±4 /iJy for the radio quiet 
sources in the three groups. We convert these measurements to the flux densities at rest frame 
5 GHz assuming a power-law radio spectrum of i v ~ jj-°- 75 (Condon 1992). This gives an 
average radio loudness (i.e. average radio to optical emission ratio (R) = (fi,,5GHz)/(fi;,44oo)) 
of (R) ~0.7, which is consistent with the typical radio loudness values of 0.1 < R < 1 found 



3 A typical redshift value of z=6 is adopted here 
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with local (Kellermann et al. 1989) and z<2 (Cirasuolo et al. 2003) radio quiet quasars. 
The derived q parameters (i.e. FIR to radio emission ratios) are all lower or at the edge 
of the range of q = 2.34 ± 0.7 for typical star forming galaxies (Yun et al. 2001), but are 
consistent with the values of q<2 found in the local Seyfert 1 galaxies (Ulvestad & Ho 2001). 
This result suggests that, as was found in the local universe and at z~2 to 4, AGN power 
dominates the radio emission in a large fraction of these radio quiet quasars, i.e., the radio 
quiet AGNs are in fact not radio silent (Barvainis et al. 2005; Carilli et al. 2001; Petric et 
al. 2006). 

We plot the mean SED (normalized by the rest frame 1450A flux density) of the radio 
quiet sources in each group in Figure 2, together with the templates of low-z optical quasars 
from Elvis et al. (1994) and Richards et al. (2006). We also plot the average optical to 
near infrared emission based on the Spitzer photometry of 13 z~6 quasars from Jiang et al. 
(2006), who show that the ratio is consistent with the local templates. The FIR-to-1450A 
ratios of the three z~6 quasar groups have a range of ~0.6 dex, and no obvious excess is 
seen between the value of the whole sample and the FIR end of the templates. However, 
due the lack of millimeter and submillimeter data, the available templates do not give a 
good constraint of the FIR-to-millimeter SED for the local optical quasars. Sanders et al. 
(1989) studied a sample of optically selected Palomar-Green quasars (PG quasars; Schmidt 
& Green 1983; Boroson & Green 1992), and their analysis suggest that the typical 100 /xm to 
1 mm spectral index a ((f^ ~ u a )) is > 2. We extrapolate the FIR emission of the templates 
with this result (an example of a = 2 is plotted in Figure 2). The extrapolation matchs the 
mean value of the z~6 MAMBO undetected quasars very well, while the average value of 
the MAMBO detection is ~2.5 times higher. 



3.2. The SEDs of objects with 250 GHz detections 

We plot the individual optical-to-radio spectral energy distribution (SEDs) for the ten 
sources that were detected by MAMBO in Figure 3. The FIR emission from the two strongest 
MAMBO detections, J1148+5251 and J0927+2001, was well measured with SHARC-II and 
SCUBA at shorter wavelengths, which show FIR emission bumps in the SEDs with dust 
temperatures of 55 K and 52 K, respectively (Robson et al. 2004; Beelen et al. 2006; Wang 
et al. 2008a). Similar excesses are also seen in the SEDs of another five sources: J0033-0125, 
J0840+5624, J1335+3533, J1425+3254, and J2054-0005, where the observed 250 GHz flux 
densities are > 0.5 dex higher than those expected from the quasar templates fit to the 
observed optical emission. 

We model the FIR emission of these sources with a optically thin graybody assuming 
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an emissivity index of (3 = 1.6 (Beelen et al. 2006). For J0927+2001 and J1148+5251, we 
directly adopt the dust temperatures fit in Wang et al. (2008a) and Beelen et al. (2006), 
respectively. For the other nine sources, we assume a dust temperature of 47K (Beelen et 
al. 2006; Wang et al. 2007). The model FIR SED was extrapolated to the radio band using 
the typical FIR to radio correlation defined by star forming galaxies, i.e. q=2.34 (Yun et 
al. 2001). A radio spectral index of -0.75 is adopted here (Condon 1992). Three out of the 
seven FIR excess sources have been detected at 1.4GHz with deep VLA observations, and 
the observed flux densities are all above the average for their 250 GHz flux and q=2.34, but 
the values are within the range in which star forming galaxies lie. 



To investigate the origin of the FIR emission, we compare the correlation between the 
FIR luminosities and the central AGN bolometric luminosities of these z~6 quasars to that of 
two samples of local quasars: a sample of forty-five local (z<0.5) PG quasars from Hao et al. 
(2005), a sample of thirty-one z<0.35 type I AGNs hosted in ULIRGs with IR luminosities 
L(8 - 1000/xm) > 1O 12 L (IR quasars; Zheng et al. 2002; Hao et al. 2005). 

Hao et al. (2005) reported that the IR quasar sample showed stronger FIR emission 
and a shallower slope in the 60 /xm-optical luminosity relationship when compared to the PG 
quasars (see also Hao et al. 2007). They attributed this result to the starburst-dominated 
FIR dust heating in the host galaxies of these IR quasars. We re-fit these FIR-AGN rela- 
tionships using the FIR and AGN bolometric (Lboi) luminosities of the two samples. The 
bolometric luminosities are taken from Hao et al. (2005), which were derived from the 
extinction-corrected quasar optical emission (see also Zheng et al. 2002). The FIR luminosi- 
ties of the IR quasars are taken from Zheng et al. (2002), while for the PG quasar sample 
we adopt the FIR luminosities and upper limits presented in Haas et al. (2003). All these 
data are corrected to our cosmology. We adopt the Expectation-Maximization method in 
the IRAF STSDAS package (Isobe et al. 1986) which can perform linear regression with 
censored data. The results are: 



Figure 4 compares the z~6 quasars to these two local quasar samples in the Lfir — Lboi 
plane. All the millimeter detected sources, and the average of the non-detections are pre- 
sented. For millimeter detections, Lfir was integrated using the modeled FIR SEDs described 



4. 



Correlations between FIR emission and the AGN 



IR quasars : logL Fm = (0.35 ± 0.08)logL bol + (7.9 ± 1.0) 
PG quasars : logL FIR = (0.66 ± 0.08)logL bol + (2.8 ± 1.0) 



(5) 
(6) 
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in Section 3.2 (see also Wang et al. 2007), and for non-detections we adopt the estimation of 
average FIR luminosity in Section 3. The AGN bolometric luminosities of these z~6 quasars 
are estimated based on optical and near-infrared observations. Jiang et al. (2006) presented 
Lboi for the Spitzer observed sources. We adopted their results (corrected to the cosmology 
we adopted in this paper) and estimated Lboi from the optical B-band luminosity for the 
other sources, Lboi = 10.4Lb (Richards et al. 2006). The B-band luminosities are estimated 
with the rest frame 1450 A magnitudes assuming a power- law spectrum of i v ~ z/~°- 5 . We 
also convert L bol to black hole accretion rate M using 

L bol = f]Mc 2 , (7) 

and an assumed efficiency r\ = 0.1. The derived parameters for the ten sources are listed in 
Table 4. 

Nine of the ten MAMBO-detected z~6 quasars stand close to the relationship derived 
from local IR quasars in the L F ir — L bo i plot (Figure 4), while the other source J0818+0722 
has a lower FIR luminosity, falling between the relationships of the IR and PG quasars. This 
may suggest that the strong millimeter detections at z~6 are the high-mass counterparts of 
the local IR selected quasars, and a similar starburst dust-heating mechanism dominates 
their FIR emission. However, these z~6 quasars are much more luminous in the optical 
than the two local quasars and the Lfir — Lboi trends derived with the two local samples are 
getting mixed at high luminosity end. Thus extending the millimeter study to the optically 
fainter quasar population at z~6 will help to determine the L FIR — L bol correlation for the 
millimeter bright z~6 quasars and make a better comparison at low and high redshifts. On 
the other hand, the average value of the MAMBO non-detections is consistent with the trend 
defined by the local PG quasars. 



5. A link between millimeter detection and the strength of Lya emission 

Omont et al. (1996b) studied a sample of z~4 optically bright quasars and found 
that the millimeter detected sources tend to have weaker UV emission lines compared to 
the average line strength of the millimeter non-detections. Bertoldi et al. (2003a) found a 
similar trend with MAMBO observations of five z~6 quasars. In this section, we investigate 
this effect with the Lya emission of full sample of z~6 quasars observed with MAMBO. 
The rest-frame Lya equivalent widths (EWL ya ) of five quasars in the SDSS deep strip were 
presented in Jiang et al. (2007a). The UV spectra of another 20 sources are available in the 
papers of Fan et al. (2000, 2001a, 2003, 2004, 2006a, 2008, in prep). We estimate EW Lya for 
18 of them. The spectra have been corrected for Galactic extinction, adopting the extinction 
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curve presented in Savage & Mathis (1979). We then fit the spectra using a method similar 
to that used in Jiang et al. (2007a), i.e. a power-law (f^ ~ z/~ - 5 ) for the UV continuum and 
gaussians for the Lya and NV line emission. The sources J1044-0125 and J1048+4637 have 
broad absorption line features and their equivalent widths are therefore unavailable. 

We plot the histograms of EWL y « for MAMBO detections (solid line) and non-detections 
(dotted line) in Figure 5. Most of the MAMBO detections distribute in the region of 
log EWL yQ (A) < 1.5 with a median value of 24 A, while most non-detections have log EWL ya (A) 
> 1.5 with a median value of 68 A which is comparable to the average value of ~70 A from 
the large optical quasar samples at lower redshifts (Schneider et al. 1991; Fan et al. 2001b). 
However, there is also one MAMBO non-detected source, J1621+5155, which shows very 
weak Lya emission with EWL ya < 5A. 

The physical origin of this curious trend is not clear. According to the spectra study 
with samples of z~4 quasars presented in Omont et al. (1996b), similar weaker emission is 
also seen in other UV emission lines, such as the Si IV/O IV] and C IV lines. Though the 
UV line emission is weak in the MAMBO detected sources, there is no evidence of strong 
dust reddening in the UV continuum. We can speculate that intrinsically different physical 
conditions associated with the broad line emission cloud, or a special dust obscuration geom- 
etry that affects only the broad line region, may explain this effect. Further tests with other 
strong broad lines, such as the C IV, and Mg II lines in the UV, the Ha, and B./3 lines in the 
optical, and the Paschen lines in near-infrared, are required to give a better understand of 
this trend. In particular, the trend may disappear when tested with the near-infrared lines 
if it is due to dust absorption. 



6. Discussion 

In this paper, we analyze the FIR properties of the current quasar sample at z~6. These 
quasars are mainly optically selected, and represent the luminous end of the quasar popu- 
lation at this early epoch. The X-ray to near-infrared SEDs of these sources are dominated 
by powerful AGN, just as in optical quasars at lower redshift (Jiang et al. 2006). We divide 
the sample into MAMBO 250 GHz detections of 10 sources and non-detections of 23 sources. 
Their FIR emission was compared to that of the low-z quasars, and with the results listed 
above we will further discuss the FIR dust heating and evolution stage of these z~6 quasars 
in this section. 
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6.1. The average FIR emission of the 250 GHz non-detections 

The average FIR emission of the MAMBO non-detections is recovered at 4<r by stacking 
the 250 GHz measurements of each source, and the mean SED is similar to the low-z quasar 
templates of Elvis et al. (1994) and Richards et al. (2006). This suggests that the emission 
from the outer part of AGN-heated dust torus may dominates the FIR SED in these sources, 
as was found in the low-z optically bright quasars (Sanders et al. 1989; Haas et al. 2003). 
However, we cannot rule out contributions from star formation in the quasar hosts. Further 
surveys of PAH and other star forming features will address the average contribution from 
star formation to FIR dust heating in these objects. The average FIR luminosity (42.5/xm - 
122.5/im, rest frame) derived from the stacking flux density is 1.2xl0 12 L . If star formation 
is not the dominating dust heating source (i.e. contributes <50%), the average star formation 
rate in the quasar hostfl should be less than 200 M yr 1 . 

6.2. The FIR emission of the strong 250 GHz detections 

Excess FIR emission is seen in seven of the MAMBO detected sources (see Figure 3). 
Their 250 GHz flux densities are much stronger (> 0.5 dex) than the values expected from 
the quasar templates when normalized in the rest-frame UV. These sources show FIR-to- 
radio luminosity ratios within the range defined by typical star forming galaxies (Section 
3.2, see also Carilli et al. 2004; Wang et al. 2007). The relation between FIR and AGN 
bolometric luminosity for these sources follows the trend defined by a sample of low-redshift 
IR luminous quasars which are hosted by ULIRGs (Section 4). The Spitzer near-infrared 
photometry by Jiang et al. (2006) detected the AGN-heated hot dust (~1000 K) in some 
of these bright millimeter sources; the optical-to-near infrared SEDs of these sources match 
the low-z templates very well. The strong FIR emission bump in these sources indicates 
an additional warm dust component, with a dust temperature of 50 to 60 K. Brightness 
temperature arguments imply that the warm dust emission region must be extended over a 
few kpc (Wang et al. 2008a). Thus either a significant fraction of the dust heating is due to 
star formation in the host galaxy, or the AGN must heat dust over large scales. 

Intense star formation activity is suggested by other observations of these strong mil- 
limeter quasars. Molecular CO emission is detected in the three strongest z~6 millimeter 
sources: J1148+5251 (Bertoldi et al. 2003b; Walter et al. 2003, 2004), J0927+2001 (Carilli 



The star formation rate is estimated with the empirical relationship from Kennicutt (1998), assuming 
a standard Salpeter initial mass function, i.e. SFR ~ 4.5Lir M^yr -1 , where Lir, is the infrared luminosity 
(8-1000 //m) in unit of 10 44 ergs _1 , and is ~ 1.5Lfir for warm dust emission. 
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et al. 2007), and J0840+5624 (Wang et al. 2008b, in prep). These CO detections indicate 
the presence of > 10 10 M of molecular gas in the quasar hosts (Walter et al. 2003; Bertoldi 
et al. 2003b; Carilli et al. 2007), which can provide the fuel required for star formation rates 
of up to > lO 3 M yr _1 . Indeed, it is hard to see how the host galaxies can avoid massive 
star formation with so much molecular gas. The measurement of several CO transitions in 
J1148+5251 suggests a gas density of about 10 4 5 cm' 3 and a kinetic temperature of ~100 
K (Bertoldi et al. 2003b), which are similar to values in starburst environments found in 
other quasars and submillimeter galaxies at lower redshifts (Carilli et al. 2002; Solomon & 
Vanden Bout 2005; Riechers et al. 2006). In addition, recent high resolution observations 
of [C II] in J1148+5251 (Walter et al. 2008, in prep.; Carilli 2008) show an emission region 
extended over many kpc, i.e., comparable to that of the CO (3-2) line. 

Further studies also suggest connections between the CO and FIR emission. The 
FIR-to-CO luminosity ratios of J1148+5251 and J0927+2001 are about 500 and 650 L 
(Kkms -1 pc 2 ) -1 , respectively, which are consistent with the values found in other high-z 
CO emitting galaxies (Solomon & Vanden Bout 2005). In Figure 6, we indicate these two 
CO detected z~6 quasars (as stars) on the FIR-CO luminosity correlation plot discussed in 
Riechers et al. (2006). The two points fall a little above but within the observed scatters 
of the FIR-CO luminosity relation of the typical star forming galaxies at both low and high 
redshifts. 

All these facts suggest that massive star formation is an important source of the FIR 
dust heating. However, it is difficult to constrain what fraction of the FIR dust heating is 
contributed by star formation in these FIR luminous quasars. According to the FIR-to-CO 
luminosity correlation in Figure 6, at least 30 — 50% of the FIR emission should originate 
from star formation, giving a corresponding star formation rate of > 1000 M Q yr _1 . 

6.3. The evolutionary stage of the z~6 quasars 

Observations from optical to near-infrared of these z~6 quasars indicate that the central 
SMBHs have masses of > 1O 9 M (Willott et al. 2003; Jiang et al. 2006). Based on the 
SMBH-bulge mass relationship (eg. Marconi & Hunt 2003), the final stellar bulge mass 
should be on order of 10 12 M for a mature system. The formation of such a massive stellar 
system may require a star formation rate of > 10 3 M yr _1 on time scales of a few hundred 
million years as was suggested by some numerical simulations (eg. Li et al. 2007a, b). 
Thus the average star formation rate of < 200 M yr _1 found in the 250 GHz non-detected 
quasars can be either the case that the major bulge building has already been finished and 
the SMBH-bulge mass relationship seen at z=0 has already been established, or that the 
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host galaxy is has not yet reached the massive star formation phase, and the growth of the 
stellar bulge is lagging behind that of the SMBH (see eg. Shields et al. 2006; Ho 2007). Deep 
and high-resolution imaging at near-IR wavelengths may provide a direct measurement of 
the stellar population in these z~6 quasar hosts. 

For the strong millimeter quasars at z~6, the FIR excess suggests a star formation rate of 
> 1000 M Q yr^ 1 from the host galaxies. It is likely that we are witnessing simultaneous build- 
ing of the stellar bulge and growth of the SMBH in these objects. The case of J1148+5251 
is the best documented: Willott et al. (2003) derived an estimate of 3 x 10 9 M for the 
black hole mass, while Walter et al. (2004) estimated a dynamical mass of ~ 5 x 1O 1O M . 
The SMBH-bulge mass ratio is thus clearly more than one order of magnitude larger than 
the local SMBH-bulge relationship (eg. Marconi & Hunt 2003), as in several other high-z 
CO-detected QSOs (Shields et al. 2006; Ho 2007). However, the amount of gas detected in 
CO observations such as in J1148+5251, a few 1O 1O M , is far short of the value required to 
reach the bulge mass given by the local black hole-bulge relationship. 

To fully understand the evolutionary stage of these strong millimeter quasars at z~6, 
further observations are required. First, CO and [C II] searches and excitation studies should 
be extended to all the z~6 quasars with strong FIR excesses to determine the host galaxy 
ISM properties. Second, high resolution (< 0.3") CO, [C II] and dust mapping should be 
made of all the CO detections. If the FIR-emitting dust and molecular gas are distributed 
similarly, the FIR emission is likely to be due to a starburst, while more compact FIR 
emission will argue for AGN dust heating. The resolved CO emission would also provide 
an estimate of the dynamical mass of the quasar hosts on kpc scales. This is a direct, and 
perhaps the only, way to test the black hole-bulge relationship at the highest redshifts. These 
observations can be done with sensitive submillimeter, millimeter and radio telescopes such 
as the PdBI, CARMA, and SMA, and the coming ALMA and EVLA. 

7. Conclusion 

We study the SEDs of 33 quasars at z~6 from FIR to radio wavelengths. We conclude 
that, when averaged with the whole sample and the 250 GHz undetected sources, no signifi- 
cant difference is seen between the mean FIR-to-radio SEDs and the templates of low-redshift 
optically selected quasars. In particular, we extrapolate the FIR emission of the templates 
with the typical quasar FIR-to-millimeter spectrum (i.e. f ~ u a , a > 2) from Sanders et 
al. (1989). This extrapolation is consistent with the average FIR emission of the 250 GHz 
undetected z~6 quasars very well (see Section 3.1 and Figure 2). This fact suggests that 
these sources have a similar AGN dominated dust heating mechanism. The average star 
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formation rate from the host galaxies of the quasars undetected by MAMBO is estimated to 
be less than a few hundred M yr _1 . 

We detect a strong FIR excess in seven of the ten z~6 quasars detected in the millimeter 
band. These FIR luminous sources are likely to be the high-mass counterparts of local IR 
luminous quasars which are in transition between the starburst phase and the mature quasar 
phase. Star formation at rates of > 10 3 M yr _1 may be the dominant mechanism to heat the 
dust that gives rise to the FIR emission. The millimeter detected sources also show weaker 
Lya emission compared to that of the non-detections, but the origin of this trend is not clear 
yet. Further observations of these FIR luminous quasars, especially spatially resolved line 
and dust studies, are required to further address questions of star formation and FIR dust 
heating in these extreme objects. 

The current sample of z~6 quasars is optically selected. There may be other quasar 
populations at z>6 which are still buried in their starburst environment and are obscured 
in the optical. Sensitive IR and millimeter facilities, such as Spitzer, ALMA, and Herschel, 
may discover such objects and provide a more complete view of galaxy and SMBH evolution 
in the early universe. 
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Table 1: Summary of the new observations 



Name 


redshift 


mi450 


Mi450 


f250GHz 


fl.4GHz 


Reference™ 










mjy 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


J000552.34-000655.8 


5.85 


20.23 


-26.47 


0.36+0.48 


40+130 


1,*,2 


J003311. 40-012524.9 


6.13 


21.78 


-25.00 


1.13+0.36 


-27+19 


3,3,* 


J020332.39+001229.3 


5.86 


20.97 


-25.73 


1.52+0.67 


195+22 


4,*,* 


J030331. 40-001912.9 


6.07 


21.33 


-25.43 


0.23+0.51 


-85+62 


4,*,* 


J035349.76+010405.4 


6.05 


20.21 


-26.55 


1.20+0.46 


17+19 


4,*,* 


J081827.40+172251.8 


6.00 


19.34 


-27.40 


1.19+0.38 


123+12 


5,*,2 


J084119. 52+290504.5 


5.96 


19.61 


-27.12 


1.00+0.43 


43+27 


8,*,* 


J084229.23+121848.2 


6.08 


19.58 


-27.18 


0.11+0.55 


-8+19 


9,*,* 


J092721. 82+200123.7 


5.77 


19.87 


-26.81 


4.98+0.75 


50+11 


5,2,* 


J104433.04-012502.2 


5.78 


19.21 


-27.46 


1.82+0.43 


-15+24 


10,*,7 


J142516. 30+325409.0 


5.85 


20.62 


-26.08 


2.27+0.51 


20+20 


14,*,* 


J142738.59+331242.0 


6.12 


20.33 


-26.44 


0.39+0.66 


1730+131 


15,*, 15 


J160253.98+422824.9 


6.07 


19.86 


-26.90 


1.41+0.54 


60+15 


1,*,2 


J162100.70+515544.8 


5.71 


19.89 


-26.77 


0.30+0.55 


-12+21 


9,*,* 


J163033.90+401209.6 


6.05 


20.64 


-26.12 


0.80+0.60 


14+15 


11,12,* 


J164121. 64+375520.5 


6.04 


21.30 


-25.45 


0.08+0.46 


-30+32 


3,3,* 


J205406.42-000514.8 


6.07 


20.67 


-26.09 


2.38+0.53 


17+23 


4,*,* 


J231546.36-002357.5 


6.12 


21.31 


-25.46 


0.28+0.60 


31+16 


4,*,* 


J232908.28-030158.8 


6.43 


21.65 


-25.20 


0.01+0.50 


14+22 


3,3,* 



Note - The detections are marked as boldface. 



"The three references are for the optical, 250 GHz, and 1.4 GHz data, respectively. The asterisks denote new data reported in 
this paper. 

References for Tabic 1 and 2 - (1) Fan ct al. 2004; (2) Wang ct al. 2007; (3) Willott et al. 2007; (4) Jiang ct al. 2008; (5) Fan 
ct al. 2006a; (6) Fan et al. 2001a; (7) Pctric ct al. 2003; (8) Goto 2006; (9) Fan et al. 2008, in prep. (10) Fan ct al. 2000; (11) 
Fan ct al 2003; (12) Bertoldi et al. 2003a; (13) Carilli ct al. 2004; (14) Cool ct al. 2006; (15) McGreer et al. 2006. 
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Table 2: Summary of previous observations 



Name 


rcdshift 


mi450 


Mi450 


f250GHz 


fl.4GHz 


Reference™ 










mjy 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


J000239. 39+255034.8 


5.80 


19.02 


-27.67 


0.20+0.88 


89+14 


1,2,2 


J083643. 85+005453.3 


5.81 


18.81 


-27.88 


-0.39+0.96 


1740+40 


6,7,7 


J084035.09+562419.9 


5.85 


20.04 


-26.66 


3.20+0.64 


12+9 


5,2,2 


J103027. 10+052455.0 


6.31 


19.66 


-27.16 


-1.15+1.13 


-3+20 


6,2,2 


J104845. 05+463718.3 


6.20 


19.25 


-27.55 


3.00+0.40 


7+13 


11,12,2 


J113717.73+354956.9 


6.01 


19.63 


-27.12 


0.10+1.13 


9+17 


5,2,2 


J114816. 64+525150.2 


6.42 


19.03 


-27.82 


5.00+0.60 


55+12 


11,12,13 


J125051. 93+313021. 9 


6.13 


19.64 


-27.14 


0.07+0.90 


37+21 


5,2,2 


J130608.26+035626.3 


6.02 


19.55 


-27.19 


-1.05+1.04 


14+21 


6,7,7 


J133550.81+353315.8 


5.95 


19.89 


-26.84 


2.34+0.50 


35+10 


5,2,2 


J141111. 29+121737.4 


5.93 


19.97 


-26.75 


1.00+0.62 


61+16 


1,2,2 


J143611. 74+500706.9 


5.83 


20.16 


-26.54 


-0.21+1.14 


6+16 


5,2,2 


J150941. 78-174926.8 


6.12 


19.82 


-26.95 


1.00+0.46 




3,3,- 


J162331. 81+311200.5 


6.25 


20.13 


-26.67 


0.17+0.80 


24+31 


1,2,2 



Note - The detections are marked as boldface. 



a The three references are for the optical, 250 GHz, and 1.4 GHz data, respectively, and the corresponding literatures are listed 
at the end of Tabic 1. 



Table 3: Mean FIR and radio emission of the z~6 quasars 

Group Numbcr a < f250GHz > < LpiR > Number 15 < fi.4GHz > < Ll.4GHz > c q 







mjy 


1O 12 L 






L Hz- 1 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Whole sample 


33 


1.26+0.10 


2.9+0.2 


32 


46+3 


0.031+0.002 


1.41+0.04 


(radio quiet) 


30 


1.29+0.10 


3.0+0.2 


29 


32+3 


0.022+0.002 


1.57+0.05 


MAMBO detections 


10 


2.73+0.06 


6.4+0.1 


10 


37+4 


0.024+0.003 


1.84+0.05 


Non-detections 


23 


0.52+0.13 


1.2+0.3 


22 


56+4 


0.038+0.003 


0.94+0.11 


(radio quiet) 


20 


0.51+0.13 


1.2+0.3 


19 


27+4 


0.018+0.003 


1.24+0.13 



a Number of sources observed at 250 GHz. 
b Number of sources observed at 1.4 GHz. 

c A radio spectral index of -0.75 (Condon 1992) is adopted here to calculate the rest frame 1.4 GHz luminosity. 
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Table 4: Luminosities of the 250 GHz detected z~6 quasars 



Name 


L bol 


M 


L FIR 


1O 1J L 


M yr" 1 


1O^L 


(1) 


(2) 


(3) 


(4) 


J003311.40-012524.9 


1.4 


9.4 


2.6+0.8 


J081827.40+172251.8 


12.8 


86.3 


2.8+0.9 


J084035. 09+562419.9 


6.5 


43.6 


7.6+1.5 


J092721. 82+200123.7 


7.4 


50.0 


12.2+2.8 


J104433.04-012502.2 


13.5 


91.1 


5.3+0.8 


J104845. 05+463718.3 


14.6 


98.4 


6.9+0.9 


J114816. 64+525150.2 


18.8 


126.8 


13.4+3.4 


J133550.81+353315.8 


7.6 


51.4 


5.5+1.2 


J142516. 30+325409.0 


3.8 


25.6 


5.4+1.2 


J205406.42-000514.8 


3.8 


25.8 


6.8+1.1 



Note - M is the black hole accretion rate derived with equation (7) in Section 4. 
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Fig. 1.— Radio 1.4 GHz field images (25'.6 x 25'.6) of J0033-0125 (a) at 6".5 x 5" resolution 
(FWHM) and J1641+3755 (b) at 6" x 4". 5 resolution. Contour levels are a geometric pro- 
gression by a factor of 2, starting at 0.120 mJybeam -1 , and the crosses denote the optical 
quasar positions. 
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Fig. 1. — Continued 
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Fig. 2. — Mean FIR and radio emission (normalized by the rest frame 1450A luminosity 
density) of the radio quiet z~6 quasars. Black squares, red circles, and green triangles 
represent the average values of the whole z~6 sample, MAMBO 250 GHz detections, and 
non-detections, respectively. We also plot the average near-infrared emission (black squares) 
derived from the Spitzer photometry of 13 sources from Jiang et al. (2006). The dashed and 
solid lines are the local quasar templates from Elvis et al. (1994) and Richards et al. (2006), 
respectively. The dotted line denotes an extrapolation of the template FIR emission (Elvis 
et al. 1994) with a power law of L u ~ v 2 
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Fig. 3. — Individual UV to radio SEDs of the 250 GHz detected z~6 quasars. Squares show 
the observational data at 1450A, 3.6 jum, 4.5 /im, 5.6 /im, 8.0 /im, 350 /im, 450 /zm, 850 jum, 
1.2 mm (250 GHz), and 20 cm (1.4 GHz), and arrows denote 3a upper limits. The 3 mm 
continuum data for the sources J0927+2001 and J2054-0005 from the PdBI CO observations 
(Carilli et al. 2007; Wang et al. 2008b, in prep.) are also included. Local quasar templates 
are plotted as in Figure-2 and normalized to rest frame 1450A. The thick solid line is a 
warm dust model (see Section 3.2) normalized to the submm data and extended to the radio 
band with the typical radio-FIR correlation of star forming galaxies, i.e. q=2.34 (Yun et al. 
2001), and the dotted lines factors of 5 excesses above and below the typical q value. 
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Fig. 3. — Continued 
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Fig. 4. — The LpiR-Lboi correlation. The 250 GHz detected z~6 quasars are plotted as 
black open squares with error bars denoting la rms. The average value of the 250 GHz 
non-detections is plotted as a gray square. The local IR and PG quasars from Hao et al. 
(2005) are plotted as circles and diamonds with arrows denoting upper limits in Lfir. The 
dashed lines represent the linear regression results for the two local quasar samples. We 
convert L bo i to the mass accretion rate (M) at the top abscissa. 
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Fig. 5. — The distribution of rest-frame equivalent widths of Lya (EWi yQ ) for MAMBO 
detected (solid line) and non-detected (dotted line) quasars at z~6. The shadowed area de- 
notes the two sources that only have upper limits for the Lya emission. The solid and dotted 
lines represent the median values for MAMBO detections and non-detections, respectively. 
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Fig. 6. — L FIR vs. CO luminosity (L' co ) for galaxies and quasars at various redshifts taken 
from Riechers et al. (2006). The dashed line represents the relationship L F ir oc Lqq 1 ' 39 , 
derived from low-z ULIRGs, starbursts, high-z submm and radio galaxies, and low-z and 
high-z CO detected quasars in Riechers et al. (2006). The stars represent the two CO 
detected quasars in the z~6 sample. 



